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^ '. ABSTRACT 

H ' We report the H^'^CO'*' (J=l-0) survey observations toward embedded clus- 



ters obtained using the Nobeyama 45in telescope, which were performed to follow 
up our previous study in the C^^O survey with a dense gas tracer. Our aim is to 
address the evolution of cluster-forming clumps. We observed the same 14 clus- 
ters in C^^O, which are located at distances from 0.3 to 2.1 kpc with 27" resolution 
^ ' (corresponding to Jeans length for most of our targets) in H^'^CO^. We detected 

O ' the 13 clumps in H^'^CO^ line emission and obtained the physical parameters 

t^^ ■ of the clumps with radii of 0.24-0.75 pc, masses of 100-1400 Mq, and velocity 

Q ' widths in FWHM of 1.5-4.0 km s^ . The mean density is ~ 3.9 x 10^ cm~^ 

and the equivalent Jeans length is ~ 0.13pc at 20 K. We classified the H^^CO^ 
clumps into three types, Type A, Type B, and Type C according to the relative 
^ ■ locations of the H^^CO"*" clumps and the clusters (see our previous study). Our 

c^ ■ classification represents an evolutionary trend of cluster-forming clumps because 

dense clumps are expected to be converted into stellar constituents, or dispersed 
by stellar activities. We found a similar but clearer trend than our previous re- 
sults for derived star formation efficiencies to increase from Type y4 to C in the 
H^'^CO^ data, and for the dense gas regions within the clumps traced by H^^CO"*' 
to be sensitive to the physical evolution of clump-cluster systems. In addition, we 
found that four out of 13 H^^CO"^ clumps which we named "DVSOs" (Distinct 
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Velocity Structure Objects) have distinct velocity gradients at the central parts 
of them, i.e., at the location of the embedded clusters. Assuming that the ve- 
locity gradients represent the rigid-like rotation of the clumps, we calculated the 
virial parameter of the H^^CO^ clumps by taking into account the contribution 
of rotation, and found that the DVSOs tend to be gravitationally unbound. In 
order to explain the above physical properties for DVSOs in a consistent way, 
we propose a clump-clump collision model as a possible mechanism for triggering 
formation of clusters. 

Subject headings: ISMxlouds — radio lines:stars — stars:formation 



1. INTRODUCTION 

Almost all stars (more than 90 % of stars within our galaxy), in pa rticular massive star s 



(> 8 Mq) within the disk of the Milky Way form as members of clusters (JLada fc Lada 20031 ). 
Thus, clusters have been long considered as important laboratories for astronomy. The 
mechanisms of cluster formation triggered by external effects have been proposed from 
both observations an d theoretical studies . For example, th e effects from H ll regions (e.g.. 



Hester fc Desch 20051) . supernovae (e.g., lYasui et al. 20061 ). and molecular outflows (e.g 



Shimajiri et al. 20081 ) toward the cluster forming regions have been well studied previously. 
In addition, Furukawa et al. (2009) argued that cluster formation can be triggered by the 
dynamical interaction of massive molecular clouds. However, evidence for triggered cluster 
formation via cloud-cloud interactions is limited to morphological evidence for most cases. 

In recent studies, the cluster-forming clumps, which are the dense regions in molecular 
clouds (size ~ 1 pc, mass ~ 100-1000 Mf^, dens ity ~ 10^~^ cm~^) are considered to be the 



parental objects of clusters ( iLada &: Lada 20031 ). Ridge et al. (2003) have carried out a 
survey of the dense clumps associated with the embedded clusters using the ^^CO(J=1-0), 
C^^O(J=1-0), and 0^*^0(7=2-1) molecular lines. They proposed a scheme for classifying 
the clumps and revealed that the peak positions of the clumps in the evolved stages are 
apart from the Herbig Ae/Be stars or IRAS sources. However, the spatial resolution of their 
observations was insufficient to resolve the components whose size is comparable to the Jeans 
length. Higuchi et al. (2009) carried out a survey of the dense clumps associated with 14 
embedded clusters in high resolution C^^O (J=l-0) observations with the Nobeyama 45m 
telescope. From the spatial relation between the clump structures and forming clusters, 
they classified the C^^O clumps into three types {Type A, i?,and C). They revealed the 
evolutionary sequence from Type AtoC with dispersing gaseous components of the clumps. 
Moreover, they found that the star formation efficiencies (SFEs) of the dense clumps have a 
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tendency to increase from Type A to C, which supports the notion that their classification 
shows the evolutionary stages of the cluster-forming clumps. 

In this paper, we present high-resolution observations with the Nobeyama 45m telescope 
toward the 14 clusters. We used H^^CO^(J=1-0) molecular line emission whose critical 
density (~ 10^ cm~^) is an order of magnitude larger than that of the C^^O line. Ikeda et al. 
(2007) found a good correlation of the spatial distribution of the dust continuum emission 
with that of the H^'^CO'^( J =1-0) line emission and identified star-forming dense cores (size 
~ 0.1 pc, mass ~ 10 M©, density ~ 10^ cm'^) with the B^^CO^{J=l-Q) line. The Hi3C0+ 
line emission traces the denser regions within the clumps associated with the embedded 
clusters. Out of our 14 targets, six are common with the sample in the study of Ridge et al. 
(2003), which are located within only Ikpc. We morphologically classify the cluster-forming 
clumps according to the spatial relation between the clump structures and the location of 
clusters. We also discuss the evolution of the dense clumps and the formation mechanism of 
the clusters. 



2. OBSERVATIONS 

We have carried out H^^CO^(J=1-0; 86.754330 GHz) observations using the 45m tele- 
scope at the Nobeyama Radio Observatory (NRO) from December 2005 to May 2008 toward 
the same targets as Higuchi et al. (2009). We used the 25-BEam Array Receiver System 
(BEARS) for the front end receiver, which had 5x5 b eams separated by 41". 1 on the plane 



of the sky (jSunada et al. 2000l : lYamaguchi et al. 20001 ). At 87 GHz, the main beam size was 



18" in FWHM and the main beam efficiency {rj) was 0.51. At the back end, we used 25 sets of 
1024 chann el auto-correlator s (AC), which had 32 MHz bandwidth and 37.8 kHz frequency 



resolution (ISorai et al. 2000l ). The frequency resolution corresponds to the velocity resolu- 



tion of 0.13 km s ^ at 87 GHz. Our ma pping observations were made using the on-the-fiy 



mapping technique (ISawada et al. 2008[ ) and covered approximately 5' x 5'-10' x 10' areas 
which corresponds to 3.2 pc x 3.2 pc at a distance of 2.1 kpc centered on the position of the 
clusters. During the observations, the double sideband system noise temperatures were in 
the range between 250 K and 350 K. We used the emission-free area near the observed sources 
as the off positions. The standard chopper wheel method was used to calibrate the intensity 
into T^, the antenna temperature corrected mainly for the atmospheric attenuation. The 
telescope pointing was checked every 1.5 hr by observing SiO maser sources near (within 
~ 45' for all target objects) the target objects. The pointing accuracy was better than 3". 
To correct the difference in intensity scales among the 25 beams of the BEARS, we used 
calibration data obtained from the observations toward the center of S140 region using the 
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SIS receiver equipped with a single sideband (SSB) filter, "SlOO", and the acousto-optical 
spectrometer (AOS) as a backend. 



We used a convolution scheme with a spheroidal function (ISawada et al. 2008[ ) to cal- 
culate the intensity at each grid point of the final map-cube (R.A., decl., and i'lsr) data 
with a spatial grid size of 9", which is half of the beamsize; the final effective resolution is 
27" and the effective integration time at each grid point is about 160 sec. As a result, the 
mean rms noise level for all of the objects were 0.15 K in TJ with a velocity resolution of 
0.13 km s-^ 



3. RESULTS 

3.1. Total integrated intensity and velocity dispersion maps of the H^^CO^ 

emission 

The left panels of Figures [THS] show the H^^CO^ total integrated intensity maps su- 
perposed on the Two Micron All Sky Survey (2MASS) images (J, if, and Kg composite 
color images). We defined the gas components traced by the H^^CO^(J=1-0) line emis- 
sion above the 5 a noise level in the individual maps as the H^'^CO^ clumps. Higuchi et 
al. (2009) detected the C^^O emission for all of the 14 objects observed in H^'^CO'*' in 
this work. For Gem 1, no H^^CO'*' line emission was detected above the 5 a noise level, 
and hence it will not be shown in our discussion hereafter. The detailed derivations of the 
physical properties are described in Section 13. 2[ The H^^CO'*' clumps show various mor- 
phologies, compact or extended, spherical or filamentary. Moreover, our high-resolution 
maps show that most of the clumps have internal structures which consist of several emis- 
sion components. We defined the internal components whose separations are comparable 
to the order of Jeans length as the sub-clumps|^ The Jeans length, Aj, was calculated as 
Aj = O.OSpc (Tk/20K) ' (?2h2/10^ cm^^) , where Tk is the kinetic temperature and nnj 
is the density of the gas components. We derive the kinetic temperature, Tk in Section 13.21 
The density of the clump, nna, was calculated as nna = (3/47r)(Mciump/ywniH2-Rciump'^), where 
^ciump is LTE mass of the clump, /imnj is mean mass of molecule, and -Rciump is the radius 
of the clump derived in Section 13.21 The definitions of peaks, cavities, and offset scales are 
given in Higuchi et al. (2009). For each target source, we defined the cluster center as the 
position of the most massive star, the peak position as position where the H^^CO'*' column 
density is locally enhanced by a factor of a few compared to the mean value of the clump 



^We defined the internal substructures as sub-clumps using the same method as Higuchi et al. (2009). 
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(H^^CO^ line emission above ~ 15 a region in all objects). We can see a variation in the 
relative positions of the stellar clusters in the distributions of the H^^CO'*' emission. For 
example, while for AFGL 5142 and S140 the cluster members are just located at the single 
peak position of the H^^CO^ emission, for Mon R2 and NGC 7129 the stellar components 
are not located at the emission peaks but at the emission skirts or valleys. The middle- 
right panels of Figures [1] to |5] show the intensity- weighted 1st moment maps of the H^^CO'*' 
emission. Some of the H^^CO'*' clumps are found to have relatively clear velocity gradients 
within the clumps, especially for S87E, S88B, AFGL 5180, AFGL 5142, Mon R2, and Ser- 
pens SVS2. There appears to be a correlation between the direction of the velocity gradients 
and the orientation of the elongated clumps. The right panels of Figures [1] to [5] show the 
intensity-weighted 2nd moment maps in H^^CO^. Most of our target cluster-forming regions 
were found to have molecular outflows in previous studies (e.g.. Wolf et al. 1990 for Mon 
R2) as shown in the right panels of Figure [T] to |3 The references for the outflow data are 
given in the captions of Figure [1] to [51 



3.2. Physical properties of the H^'^CO^ clumps 

In this section we describe the derivation of the physical properties of the H^^CO'*' 
clumps. In Table [H we summarize the physical properties derived using the following meth- 
ods. In order to define the clump radius, -Rciump, we firstly calculated the observed radius, 
-Robs- We directly measured the projected area of the H^^CO^ clump on the sky plane A, 
which is the emission area defined by the 5 a contour (mean intensity is ~ 0.8 K). The 
observed radius was derived as i?obs=(^/7r)^^^- -Rciump "was estimated with correcting for 
broadening effect by the observing beam as 



-Rclump = [-Robs ^ (^FWHm/2) ] ' , (1) 

where 6'fwhm of 27" is the effective beam size in FWHM. 

The velocity width of the clump (one-dimensional velocity widths in FWHM), AV^iump 
was derived as follows. We derived the H^^CO"^(J=1-0) line widths (AVobs in FWHM) by 
fitting the H^^CO^(J=1-0) spectra averaged over the clumps with a Gaussian function. We 
estimated intrinsic velocity widths of the clumps, AKimnp, by correcting for the line broad- 
ening by profile synthesis in auto-correlators, AV^pec, as AV^iump =(^Kbs ~ ^Kpcc) • We 
found that the velocity widths, AV^iump range from 1.5 km s^ to 4.0 km s~ , which indicates 
that non-thermal contribution is dominant; the thermal velocity widths (one-dimensional 
velocity widths in FWHM) of Hi3co+ molecules are 0.08-0. 14 km s"^ with Tk = 10-30 K. 
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The H2 column density of the clumps, A^(H2) was calculated as, 

Xui3r'n+ \ I To^\ I 4.16 



^(«^)-i-°xi''"(5lfTFiT)"(x)-p( 



Te^/K 



x(^)-7 V^V^^) [cm-], (2) 

VO.Siy Vl-exp(-r)y VKkms~V 

where Xjji3co+ i^ the fractional abundance of H^'^CO^ relative to H2, Tgx is the excitation 
temperature of the transition, r is the optical depth, and / T^dv is the total integrated inten - 



sity of the H^^CO'*' line emission. For Xjji3co+) "^^ used 5.6 x 10^^^ ( JAoyama et al. 200ll ). 
Here, we applied the optically thin condition (r < 1) in H^^CO"*" line emission because we 
found that the H^^CO'*' line emission peak is weak er than the peak int ensity of C^^O line 



emission which is considered to be optically thin ( iHiguchi et al. 20091 ) . In order to esti- 
mate the excitation temperature of H^^CO^, Tgx, we used the NH3 data from the NRO 
data archive, which have already been presented in Higuchi et al. (2009). Under the local 
thermodynamic equilibrium (LTE) condition, we consider that the excitation temperature 
is comparable to the kinetic temperature (Tgx ~ Tk). We applied the relation between the 
kinetic temperature and the rotation temperature, Tk ~ Trot as proposed by Danby et al. 
(1988). The rotation temperatures, Trot(2,2;l,l) were derived from the NH3 data using the 
same method presented by Ho & Townes (1983). 



The total mass under the LTE condition, Mdump is calculated as, 

^clump = 15 



>^[-^VM^Pk-A[7i^] (t/t^^I Iau (3) 



5.6x10-11; V2100pcy V 50" ) 

Kj ^ V^ex/K/ V0.51/ VKkms" 

where D is the distance to the object and -Rciump is the radius of the clump. Higuchi et al. 
(2009) calculated the virial mass of the clumps under the condition that the internal energy 
of the clumps is contributed by the thermal and turbulent motions of cloud constituent. 
However, this condition is likely invalid for the Hi^CO^ clumps because some of those have 
distinct velocity gradients. Instead, we discuss the energy balance of the Hi^CO^ clumps by 
taking into account the contribution of rotation (see Section | 



3.3. Comparison between the physical properties of the Qi^O clumps and the 

H13C0+ clumps 

In this section, we focus on the spatial relations, and the comparison of the physical 
parameters between the Qi^O and Hi^CO^ clumps. The middle-left panels of Figure [1] to 
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|5]show the total integrated intensity maps of the C^^O(J=1-0) emission (contours; Higuchi 
et al. 2009) superposed on the H^'^CO^ emission. By comparing the spatial distribution 
between C^^O and H^^CO'*' line emission in each object, we found that H^'^CO'*' emission is 
concentrated where the C^^O line emission is strong and the H^^CO"*" clumps are found to 
have compact and simpler shapes than the C^^O clumps, which suggests that the H^'^CO^ line 
emission traces dense region within the C^^O clumps. For example, there are very compact 
H13C0+ clumps within the extended C^^O clumps in S88B. Mon R2 and NGC 7129 which 
have complex structures in C^^O emission have compact multiple emission components in 
H^^CO'*'. For most of the observed objects, the positions of the emission peak of H^^CO^ 
agree with those of C^^O. Some of them, S87E and Serpens SVS2 have emission peaks of 
H^^CO^ which are in disagreement with those of 0^*^© by an order of Jeans length. 

By comparing the LTE mass between the 0^*^© and the H^^CO^ clumps, we found that 
the ratio of Mciump(H^'^CO^) to Mciump(C^^O) is below unity in 90% of our targets. Only 
for AFGL 5180, the ratio of Mciump(H^3C0+) to Mciuinp(C^*^0) becomes unity. Comparison 
of the velocity widths between C^^O clumps and H^'^CO'*' clumps suggests that the ratio of 
Ar(Hi3C0+) to AV"(CiS0) is below unity in half of our targets. The rest (50%) of our 
targets with relatively massive stars tend to have larger ratios than unity. They have the 
large AV{}1^'^C0^) (> 2 km s'^) values for S87E, S88B, Gem 4, AFGL5180, AFGL5142, and 
S140. Two thirds (except for Gem 4, S140) of our target clumps with large Ay(H^'^CO^) 
have clear velocity gradients in the 1st moment maps. The detailed parameters are listed in 
Table H 



4. DISCUSSION 

We observed the H^'^CO^ emission toward 14 embedded clusters in order to address 
the formation and evolution of clusters. In this section, we discuss the evolutionary status 
and the formation mechanism of clusters from the observed physical properties in terms of 
the following aspects: (1) evolution of dense gas and star formation efficiency, (2) kinematic 
structure of dense clumps, and (3) clump-clump collision as a possible triggering of cluster 
formation. 



4.1. Evolution of dense gas in cluster- forming regions 

4.1.1. Morphological classification of H^^CO^ clumps 

In this section, we focus on the morphological features of H^^CO^ clumps and the 
relation of the spatial distribution between H^^CO^ clumps and stellar clusters. Our obser- 
vations revealed that the H^^CO^ clumps show various morphologies, compact (e.g., AFGL 
5142) or extended structure (e.g.,Mon R2), cometary (e.g.,S140) or filamentary (e.g.. Gem 
4, Serpens SVS2). Most of them have single or double components within the clump. We 
categorize the H^^CO^ peak positions into two cases, "Peak" and "Off-peak", using the 
separation between the H^^CO'*' peak positions and the cluster centers ; if the separation is 
smaller than the Jean length estimated for each target, we call them the "Peak" or instead of 
the "Off-peak" . The number of identified sub-clumps is dependent on the spatial resolution 
of the observations and the distance of the target objects. Note that for a distant object, 
two or more sub-clumps might be identified as one component of sub-clump. Nevertheless, 
for S87E at the distance of 2.1 kpc, the separation between the sub-clumps corresponds to 
~ 0.3 pc which is twice of Jeans length. Hence, the above discussion still holds with our 
definition and identification of the sub-clumps. Additionally, the H^^CO^ maps show a vari- 
ation of brightness contrast. Under the optically thin condition for the H^'^CO^ line, we 
can assume that the brightness contrast is equivalent to the column density contrast, which 
reflects the physical state of the dense regions within the clumps. There is a variation in 
the relation of the spatial distribution between H^'^CO^ peaks and cluster centers. Extreme 
cases are found in two sample groups; a stellar cluster located just at a single peak in the 
distribution of the H^^CO^ line emission (e.g., AFGL 5142, S140), or located at a cavity of 
the emission line distribution within the clump which has multiple components around the 
cavity (e.g.,Mon R2). As an intermediate case between them, a cluster is associated with 
one of the H^^CO^ peaks (e.g., AFGL 490) or not associated with H^^CO^ emission peak, 
but located within the clumps (e.g.,S87E). There are two distinct clusters in Gem 4, and 
one of them is located at the emission peak of H^^CO'*'. 

Higuchi et al. (2009) proposed a morphological classification using the C^^O emission 
line for the same targets. We applied the same classification scheme to the H^^CO^ clumps 
and sort the clumps into three types to simplify discussion and to clarify a qualitative trend. 
We define Type A to be clumps in which the clusters are just associated with a single peak 
of H^^CO^ emission distribution. These clumps not only have a single peak but a higher 
brightness contrast in the distribution of the H^^CO'*' emission. We define Type C to be 
clumps in which the clusters are located apart from the main H^^CO^ emission, or at a 
cavity-like H^^CO'*' emission hole. The clumps of this type have lower brightness contrast 
than Type A. The rest are Type B clumps in which the cluster is associated with one of the 
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peaks of the H^'^CO'*' emission distribution, or not associated with H^^CO'*' emission peak 
but located within the clumps. The clumps of this type have relatively higher brightness 
contrasts than Type C. We classified AFGL 5142 and S140 into Type A. In contrast, we 
classified S88B, S235AB, Mon R2, NGC 7129, and Serpens SVS2 into Type C. In addition, 
we classified S87E, Gem 4, AFGL 5180, GGD12-15, BD+40°4124, and AFGL 490 into 
Type B. The detailed parameters of individual sources are listed in Table [2l 

We found that the above stages of the H^^CO'*' clumps are consistent with those of 
the C^^O clumps except for AFGL 5180, and Gem 1. This result shows that the dense 
regions of the clumps are extracted by the H^^CO'*' line emission more than those traced 
by the C^*^© line emission. As in Higuchi et al. (2009), we confirmed that the evolutionary 
sequence from Type ^ to C progresses with dispersing the H^^CO'*' clumps. For AFGL 
5180, the classes are different depending on which line is used; Type C (C^®0) and Type B 
(H^^CO'*'). The reason for classifying this object as Type C is mainly because of the multiple 
peaks of the C^^O clump and the cluster located between them. On the other hand, the 
H^^CO^ clump has a single peak, thus, it is classified into Type B, which is robust for a dense 
region of cluster-forming clump because of the consistency of dis tribution in the H^^CO'*' 



emission with that in dust continuum emission (JKlein et al. 20051 ). We suggested that the 
classes which represent the evolutionary stages of cluster-forming clumps are more suitable 
to be determined from the H^^CO"'' data than from the C^^O data. The clear trend in the 
relation of classes and SFEs from the H^^CO'*' data shown in Section [4.1.21 also agrees with 
this suggestion. For Gem 1, we did not detect H^'^CO'*' emission although we detected and 
identified the C^^O clump, which was classified as Type B. From the map of dust continuum 



emission (JKlein et al. 20051 ). it is understood that the emission component of dense region is 
too weak and compact to be detected significantly even with our observations. We estimated 
an upper limit to the dense gas fraction defined by Mciump(H^'^CO''')/Mciuinp(C^^O) with the 
5 a noise level of the H^^CO'*' map, and found it to be 0.01, an order of magnitude smaller 
value than those of other objects. Thus, we consider that Gem 1 is an exceptional case, in 
which the dense gas might be regulated to form or dispersed effectively. A schematic picture 
of these three types using the spatial distributions of C^^O and H^'^CO^ is shown in Figure 
El 

4- 1-2. Star formation efficiency variation along the clump-cluster evolution 

We proposed the evolutionary scenario of dense clumps associated with embedded clus- 
ters in Higuchi et al. (2009). Since a large amount of gas is dispersed by the cluster members 
during cluster evolution, stellar mass becomes more dominant in the clump-cluster system. 



-10- 



We found a tendency along the type sequence A- C for the clump structures to become com- 
plex having some sub-clumps and local peaks. The typical separation among the sub-clumps 
is the order of the Jeans length, which supports the following interpretation. In molecular 
clumps, the column density structure is more or less uniform initially. Eventually, frag- 
ment ation occurs with so me mechanisms in the clumps and initiates star cluster formation 



[e.g. 



Bonnell et al. 20031 ). After fragmentation, the local condensations grow to become 
distinct peaks where stellar clusters will eventually be born. During the cluster formation, 
the cores within the clumps form stars in a wide mass range, and dense gas is dispersed 
by their stellar activities. Once clusters are formed, young stars, particularly the most 
massive stars, disperse surrounding ga s with evolution making cavities around them (e.g. 



Fuente et al. 19981 : iFuente et al. 20021 ). The structures of the evolved clumps are expected 



to be less condensed than those of the younger clumps. For more evolved clumps, continuous 
dispersal of the dense gas around the stellar cluster makes a disagreement between the stellar 
positions and the peaks of gas condensations in which clusters are located at the cavity-like 
structure of the clumps. 

From our proposed scenario, the star formation efficiency (SFE = Mduster/ (^cluster + 
^ciump) where Mduster is the stellar masses of clusters, which are listed in Lada & Lada (2003) 
and Mciump is the masses of clumps) is expected to increase from Type A to Type C. Only 
for S140 and BD+40°4124, the total masses of the stellar members are estimated using the 
^star-^ciuster relation that was derived using the data in Lada & Lada (2003). We calculated 
the SFEs for our targets and listed the results in Table [H We found that the SFEs for the 
H^^CO"*" clumps range from 3 % to 35 % with a mean of 20 %, which is larger than the values 
previously derived for low-mass star forming regions (e.g., Onishi et al. 1998; Tachihara et al. 
2002 for the results of C^^O line emission). Figure [7] shows a histogram of the SFEs of each 
type and indicates that Type A objects have smaller SFEs (< 8%) than Type B and Type C 
objects (> 10%). We found on a similar trend in C^^O data; the larger SFEs are derived 
along the classes from Type A to C. The velocity gradients for DVSOs appear clearer in 
H^^CO"*" than in C^^O. In previous studies, the c orrelation between changing of morphologies 



and SFE along evolut ion was suggested (e.g., iLada &: Lada 20031 : [Thompson et al. 2006 



Gutermuth et al. 20091 ). Higuchi et al. (2009) demonstrated a modest correlation between 
morphological variation and SFEs. In this paper, we confirmed that the correlation is clearer 
with H^3C0+ than C^^O (see Figure 8 (b) in Higuchi et al. (2009)). This result indicates that 
H^'^CO^ traces the dense regions in the C^^O clumps presented in Higuchi et al. (2009) and 
physical conditions of dense regions are sensitive and affected by the clump-cluster evolution. 
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4.2. Kinematic structures of the clumps 

We found that half of our objects (S87E, S88B, AFGL 5180, AFGL 5142, Mon R2, 
and Serpens SVS2) have clear velocity gradients in the 1st moment maps as shown in the 
middle-right panels of Figures [1] to [51 The others have random velocity structures in the 1st 
moment maps. For the C^^O clumps of the same targets presented in Higuchi et al. (2009), 
none of them have clear velocity gradient in the 1st moment maps. The C^^O line emission 
is considered to trace not only the dense region but also surrounding gas of the dense region, 
and to barely clarify the velocity field in the dense region. We suggest that the velocity 
gradient is one of the most important features of the dense region in the cluster-forming 
clumps for the understanding of the evolution of clumps that form stellar clusters. 

In order to investigate the velocity structures of the H^^CO^ clumps in detail, we made 
position-velocity (P-V) diagrams for six objects which have relatively clear velocity gradients 
in the 1st moment maps. The axis of the slice is parallel to the direction of velocity gradients 
(see Figure [51 to [TU]) . There are the bulk linear gradients which include the multiple compo- 
nents in P-V diagram for the four objects out of six. The velocity gradient in the H^'^CO"'' 
clumps were derived for all 13 objects with a linear fitting to the emission in the P-V dia- 
grams, which are clipped at the 3 a noise levels in the channel maps. The derived values of 
the velocity gradients, Ugrad, are listed in Table O The gradients range from 0.5 kms~^ pc~^ 
to 4.3 kms~^pc^^. The gradients for the four objects are as large as 2.2 kms^^pc^^ to 4.3 
kms'^^pc^^. In addition. Figure [TT] shows the spatial distribution of the blue-shifted and 
red-shifted components with the clusters, and all of the clusters are found to be located in 
the overlapped areas between blue-shifted and red-shifted components; a typical example 
is S87E, showing two spatially separated clumps at both sides of the cluster center. The 
velocity ranges of the blue-shifted and red-shifted components are described in the caption 
of Figure [111 The spatial distributions of the blue-shifted and red-shifted components imply 
the molecular outfiows. The molecular outfiows in our targets have already been detected 
in previous studies. We found the following features. (1) There are no correlations between 
the orientation of outfiows and those of the velocity gradients. (2) These outfiows have been 
detected in the optically thick tracers (e.g., CO lines), in contrast we expected that the line 
emission is optically thin. (3) In the spectrum data, there are no wing components. From 
the above observational evidence, we conclude that the velocity structures are not made by 
outfiow activities. There is a feature that the four objects (S87E, S88B, AFGL 5180, and 
AFGL 5142) with the above velocity structures contain relatively massive stars (> 10 Mq) 
within the clusters. Here, we named "DVSOs" (Distinct Velocity Structure Objects) after 
the objects with distinct velocity structure for the above objects. These interesting kinematic 
structures are discussed in Section 14.31 in terms of the clump-clump collision model. 
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We examined the energy balance of the clumps using the equilibrium virial theorem 
(e.g., Goodman et al. 1993). The gravitational potential energy is estimated from Q = 
— OvirGM^j^ /i?ciump5 where Ovir is the coefficient Ovir = (3— p)/(5— 2p ) with the index of den - 



sity proffie, p for the correction for derivation from constant density (IWilliams et al. 19941 ). 
In the assumption of Goodman et al. (1993), they adopted p = 2 using the isothermal 
collapsing core models of the radial profile of the core. The center position of radial proffie 
is defined as the position of star within the core. In contrast, the clump has inner compo- 
nents (e.g., sub-clumps) so that the density profile is not as simple as that of the core. It 
is difficult for us to define the center position of radial profile toward the cluster-forming 
clump. In this paper, we assume the uniform density sphere to describe the inner density 
structure of the clumps. Assuming that the linear velocity gradient in each clump represents 
rotational motion, we estimated the rotational energies from E^ot = |«rot^ciump<^^-Rciump; 
where ttrot = 2/(5 + 2p), oj is the angular velocity. We estimated the turbulent energies 
as, fi'turb = |M:iump('^^urb/81n2), and the turbulent velocity widths AV^u^b were derived 
by subtracting AVthcrm and the velocity gradients from the observed velocity widths. By 
dividing the virial equation by \Q\, we calculated /3vir = 2(/3therm + Aurb + /3rot) — 1, /3vir is 
the total energy of system divided by \Q\, and Aherm, Aurb, and (3rot are the virial ratios of 
the thermal, turbulent, and rotational energies, respectively. Especially, it is noteworthy to 
compare the virial ratios of the rotational energies, which are listed in Table [31 

From the observational data, we find /3vir to be 0.05 - 0.79 (average ~ 0.25). For 
DVSOs, the /3rot values range from 0.33 to 0.79, which are comparable to the result obtained 
by Takahashi et al. (2006), who found a core with a large velocity gradient in the OMC-3 
region. The mean /^rot for the DVSOs is a factor of 3.7 higher than the mean of the non-DVSO 
sources, not a factor of 10. We also estimated /3vir and the derived values are listed in Table 
[HI The parameter /9vir is an indicator of gravitational stability; /3vir=0 for gravitationally 
bound state, /3vir<0 for contracting, and /3vir>0 for being unbound (expanding). We found 
the obtained /3vir values range from —0.44 to 1.9 (average ~ 0.39). The objects with large 
/Srot values tend to have positive /3vir values, which result from the contribution of rotational 
energy to a gravitationally unbound condition, even if we take account of uncertainty of 
/3vir. Fifty percent of our targets have (3^i^ that are less than or close to zero, which implies 
that they are gravitationally bound. For NGC 7129 and Serpens, which are classified into 
Type C, the derived /^rot values are positive. This result suggests that the gas dispersal of 
the clump leads the gravitationally unbound conditions. 
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4.3. Clump-clump collision as a possible triggering of cluster formation 

As discussed in the previous section, we found the following kinematic signatures in 
four out of the 13 H^^CO'*' clumps (DVSOs); (1) there exist distinct velocity gradients, (2) 
they have multiple components with different velocities on the rigid-like rotation, and (3) 
the clusters are located at the overlapped areas of the blue-shifted and red-shifted compo- 
nents in Figure [TTl Thus, we regard the above kinematic condition as being one of the key 
structures in cluster formation, especially for the clusters including massive stars. We note 
the characteristics of DVSOs as follows. The velocity gradients are seen in the objects for all 
three morphological types. Although whether the velocity gradient of clumps with rotational 
motion is observable depends on the orientation of the rotational axis to the line-of-sight, 
the uniformly distributed detection in the Types of objects seems to indicate less correlation 
of these velocity structures with the three morphological types. We assume the velocity gra- 
dients as the rotational motions of the dense gas in the clumps, the contribution is a large 
portion of the internal energy within the clumps. As a result, the clumps tend to be in a 
gravitationally unbound state. Such cluster formation occurring in a gravitational unbound 
dense clump is somewhat puzzling. These observational results imply that these velocity 
structures are related to the initial conditions of cluster formation in the dense clumps, and 
must be preserved during the ages (~ 1-3 Myr) of our targets. We propose a clump-clump 
collision model as a possible mechanism for triggering formation of clusters to reproduce our 
observational results on the H^^CO^ clumps. 

One possibility that can solve this puzzle is a clump-clump collision. We interpret 
the obtained data as the post-interacted (collision) stages of the clumps. Clump-clump 
collision can be responsible for triggering cluster formation through gas compression at the 
overlapped areas of the collision; shock heating at the overlapped areas can increase Jeans 
mass (oc T^^"^/ y/p) and lead to form high-mass stars. The offset collision which is more 
common than head-on is also responsible for producing the angular momentum in merged 
clumps. Such a merged clump system is not necessary to be gravitationally bound for 
inducing cluster formation, although the collision velocity should not be too high to destroy 
the clump system. If the clump-clump collision induces the angular momentum of the 
merged clump system, the dynamical timescale of the rotation. Trot is estimated as ri.ot=27r/a;, 
where u is the angular velocity of the clump. For the typical angular velocity of DVSOs, 
a;~2xl0~^^ s~^, the timescale, r^ot is calculated to be ~ IxlO^yr, which is comparable to 
the free-fall time. If the rotational motion survives for the same period as the dynamical 
time scale, then the grad ient can be observed for Type C clumps, which have ages of ~ 1-3 



Myr (JLada fc Lada 20031 ). The remnant of the collided clumps would also survive in the 
merged clump system and would be observed as multiple velocity components in H^^CO^. 
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Recently, dynamical interaction of massive clouds to induce cluster formation has been 
proposed in observational and theoretical studies. In the previous results, physical features 
similar to our results have been found. Furukawa et al. (2009) proposed a scenario of cloud- 
cloud collision toward the cluster of Westerlund 2. They found that the Westerlund 2 is 
located between the two C0(J=2-1) clouds, although their sizes and masses are an order 
of magnitude larger than those of our observed clumps. They found that the CO clouds 
are not gravitationally bound and expanding outward from the clusters, and proposed that 
the dynamical interaction between the clouds which promoted the gravitational instabilities 
form the Westerlund 2 cluster. For S87E, which is one of our observed targets, Xue & Wu 
(2008) suggested from the distributions in the multiple wavelength maps that the cluster 
formed by the clump-clump collision. Their results show that two clumps are approaching 
and the cluster is located at the overlapped area between the two clumps. In addition, 
they suggest that the distribution of the compact H ll region supports the above scenario. 
While such observational evidence for cloud-cloud collision exists, the previous studies were 
mainly based on the spatial distribution between clouds and clusters in most cases. In the 
theoretical perspective, Bhattal et al. (1998) simulated the phenomena of collisions of two 
clumps (size ~1 pc, mass ~1OOM0, collision speed ~lkms~^) in a wide range of impact 
parameters. Their results show rigid-like rotations of the gas motion induced by clump- 
clump collision, for isothermal spheres, in the collision with a small impact parameter (the 
impact parameter is 0.3). In their simulation, the clump-clump collisions produce multiple 
stars from B6 to BO stars. The spatial and kinematic properties in our results are consistent 
with those of Bhattal et al. (1998). Their calculated /Srot sue 0.1-0.2 during the interaction. 
We also estimated /Srot values of our observed objects adopting p = 2, and the derived /Sfot 
values are consistent with the results of Bhattal et al. (1998). These consistencies indicate 
that cluster formation triggered by clump-clump collision can reasonably explain the four 
H^'^CO^ clumps with velocity gradients, forming massive cluster members. 

Another possibility is turbulent fragmentation. The kinematic signatures are pro- 
duced even in the early phase of clusters by the internal activities of the cluster mem- 
bers (e.g., outflows from newly born stars) or the external effects (e.g., nearby H ll re- 
gions, supernovae etc.), which drive the turbulence energy toward dense clumps and en- 
hance a turbulent fragmentation. This mechanism may induce multiple cloud compo- 
nents with different velocities. H owever, our targets are not contained in the catalogs of 



galactic supernovae (IGreen 20091) and bubbles, which are detected in 5'iPzfeer/GLIMP SE 



survey (IChurchwell et al. 2006r). by contrast the Cassiopeia A (e.g., iMaeda et al. 20091 ) or 



W3/W4/W5 regions (e.g.. lAllen et al. 20051 : iNiwa et al. 20091 ). There are no spatial correla 



tions between the areas of large velocity dispersion and the outflow lobes (see the right panel 
of Figure [1] to [5]) . It takes a long time to produce large velocity gradients and these effects 
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are hardly detectable unless they appear in evolved phases of cluster formation including 
an energetic internal source. This is because the efficiency to convert the outflow energy 
into the turbulence energy of the dense clumps is generally expected to be low (a few % or 
less; Dyson & Williams 1997). Therefore, it is difficult to explain that the relatively young 
clumps categorized in Type A (e.g., AFGL 5142) have the distinct velocity gradient. 

The cloud-cloud (clump-clump) collision mechanism is previously considered as an ef- 
ficient mechanism to trigger star formation, but it also leads th e destruction of the cluster- 



forming systems rather than gravitational collapse ( iGilden 19841 ) . The fraction of star forma- 
tion triggered by cloud-cloud collision may be low in our Galaxy ( lElmegreen 19981 ). However, 
our results show that there are many more possibilities, especially to form massive stars by 
the clump-clump collision mechanism, which might occur inside the giant molecular clouds. 
In previous studies, no clear velocity gradients have been detected and the velocity structure 
has not been discussed for the cluster-forming clumps. Our results suggest that it is impor- 
tant to observe more dense regions in the cluster-forming clumps and investigate the velocity 
structures to understand the formation mechanism of clusters. In addition, to understand 
the initial conditions of the cluster formation, we would need to focus on the objects which 
are considered to be in younger phases of cluster formation (e.g.. High mass protostellar 
objects, and Infrared dark clouds). 



5. CONCLUSION 

We have carried out a survey toward the 14 embedded clusters in the H^^CO^(J=1-0) 
line with the Nobeyama 45m radio telescope. The purpose is to determine the evolutionary 
stages of the cluster-forming clumps and to understand the formation mechanism of massive 
clusters from the obtained physical properties. Our results and conclusions are summarized 
as follows. 

1. We made the H^^CO"*" maps with a size of ~ 5' x 5'-10' x 10' of the 14 nearby {D < 
2.1 kpc) cluster-forming regions. As a result, 13 clusters out of our targets are found 
to be associated with the H^'^CO^ clumps, whose radii, masses, and velocity widths 
are 0.24-0.75 pc, 100-1400 M©, and 1.5-4.0 km s-\ respectively. 

2. We classified the H^'^CO^ clumps into three types ( Type A, Type B, and Type C) based 
on the spatial relation between the structures of H^^CO'*' clumps and the locations of 
clusters using the same method in Higuchi et al. (2009). The evolutionary stages 
determined from the H^^CO''' data are consistent with those from the C^®0 data, 
which implies that our classification scheme is robust. 
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3. The tendency of SFEs estimated from the H^^CO^ data to be larger along the types 
from ^4 to C is similar to that from the C^^O data. The discrimination of SFE distribu- 
tion among the evolutionary types is more conspicuous in H^^CO'*' data than in C^^O 
data. This result indicates that H^^CO'*' traces the dense regions in the C^^O clumps 
and physical conditions of dense regions are sensitive and affected by the clump-cluster 
evolution. 

4. Four out of 13 H^^CO'*' clumps, we named DVSOs have the distinct velocity gradients 
within the 1st moment maps. We found that the embedded clusters are located at 
the center of the velocity fields of dense clumps. We estimated the virial ratios taking 
into account the contribution of rotation, and found that DVSOs are likely to be 
gravitationally unbound conditions. 

5. In order to explain the physical conditions of DVSOs, we propose a clump-clump 
collision model as a possible mechanism for triggering cluster formation. 
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Fig. 1. — Total integrated intensity maps of the H^'^CO^(J=1-0) emission (contours) super- 
posed on the JHKs composite color images in log scale from 2MASS (left), total integrated 
intensity maps of the C^^O(J=1-0) emission (contours; Higuchi et al. 2009) superposed 
on the H^^CO^(J=1-0) emission (middle- left), H^'^CO^(J=1-0) emission (contours) super- 
posed on the 1st moment maps (middle-right), and velocity dispersion maps (right) for S87E 
(a), S88B (b), and Gem 4 (c). The contours with the intervals of the 5 a levels start from 
the 5 a levels, where the 1 a noise levels are 0.09 K km s~^, 0.08 K km s~^, and 0.05 K km 
s~^ in Tl for S87E, S88B, and Gem 4, respectively. The outlines of the blue-shifted and red- 
shifted outflows are plotted on the velocity dispersion maps. The filled circle at the bottom 
left corner in each panel shows the effective resolution in FWHM of 27". The references of 
outflows are Barsony (1989), Phillips & Mampaso (1991) and Snell et al. (1988) for S87E, 
S88B and Gem 4. The stellar marks indicate the position of cluster center. The pink arrows 
show 1 pc scale of individual regions. 
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Fig. 2.— Same as Figure Efor AFGL 5180 (a), GGD12-15 (b), and S235AB (c). The 1 a 
noise levels are 0.06 K km s'\ 0.07 K km s-\ and 0.04 K km s"^ in T^ for AFGL 5180, 
GGD12-15, and, S235AB respectively. The references of outflows are Little et al. (1990) 
and Nakano & Yoshida (1986) for GGD12-15 and S235AB. 
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Fig. 3.— Same as Figure [D for AFGL 5142 (a), Mon R2 (b) and S140 (c 

levels are 0.07 K km s~^, 0.07 K km s^^, and 0.12 K km s 

R2, and S140, respectively. The references of outflows are Hunter et al. (1995), Wolf et al. 

(1990) and Hayashi et al. (1987) for AFGL 5142, Mon R2 and S140. 
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Fig. 4.— Same as Figured for AFGL 490 (a), BD+40°4124 (b) and NGC 7129 (c). The 1 
a noise levels are 0.06 K km s-\ 0.08 K km s-^ and 0.07 K km s"^ in Tl for AFGL 490, 
BD+40°4124, and NGC 7129, respectively. The references of outflows are Mitchell et al. 
(1995), Palla et al. (1995) and Eiroa et al. (1998) for AFGL 490, BD+40°4124 and NGC 
7129. 
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Fig. 5. — Same as Figured] for Serpens SVS2 (a). The 1 a noise levels are 0.06 K km s ^ in 
T^. The reference of outflows is Davis et al. (1999) for Serpens. 
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TypeA 

(AFGL5142, S140) 

" single peak 

' Brightness contrasts in C^^O 

and H13C0* are high 

' Cluster is associated with the peak 

within Jeans length 




Cluster 



H13C0* clump 



Ci«0 clump 
"Peakof H'3CO+ clump 

Peak of CiaO clump 



Type B 

($87E, AFGL5180, Gem4, 
GGD12~15, BD+40"4124, 
AFGL 490) 

* Double or triple peaks in Gi*0, single or double peaks in H^^CO' 

* Brightness contrasts In C''^0 and H'^CO^ are relatively high 

* Cluster is associated with one of the C^'O peaks, 
whose separation is larger than Jeans length 





Type C 

(S88B, S235AB, 

IVIon R2, NGC 7129, Serpens) 

* Triple or multiple peaks 

* Brightness contrasts in 0^80 and H'^CO^ are low 

* Cluster tends to be apart from the peaks, 
whose separation is larger than Jeans length 

Fig. 6. — Schematic picture of our morphological classification of cluster-forming clumps. 
Type A clumps in which the clusters are just associated with a single peak of C^^O and 
H^'^CO'^ emission distribution. Type B clumps in which cluster is associated with one of 
the peaks of 0^*^© emission distribution, and associated with one of the H^^CO'*' emission 
peaks. Some of them are not associated with H^'^CO"^ emission peak, but located within 
the H^^CO^ clump. Type C clumps in which the clusters are located at a cavity-like C^^^O 
and H^^CO"*" emissio n hole. Individual sources are classified into the same stages in C^^O 
( jHiguchi et al. 20091 ) . The classifications are interpreted as the evolution of cluster- forming 
clumps from Type ^ to C. 
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Fig. 7. — Histogram of the SFEs of H^^CO^ clumps. The pink color shows Type A object, 
the yellow color shows Type B object, and the blue color shows Type C object. 
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Fig. 8. — Total integrated intensity maps of the H^'^CO^( J=l-0) emission (contours) su- 
perposed on the JHKs composite color images in log scale from 2MASS for S87E (a) and 
S88B (b), which are the same as the left panels of Figure [1] (left). The right panels are 
the position- velocity (P-V) diagrams of H^^CO'^(J=1-0) emission. The contours in the P-V 
diagrams with the intervals of the 3 a levels start from the 3 a levels, where the 1 a noise 
levels are 0.06 K, and 0.06 K in T^ for S87E and S88B. The pink arrows show 1 pc scale of 
individual regions. 
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Fig. 9.— Same as Figure E] for AFGL 5180 (a) and AFGL 5142 (b). The contours in the 
P-V diagrams with the intervals of the 3 a levels start from the 3 a levels, where the 1 a 
noise levels are 0.06 K, and 0.08 K in T^ for AFGL 5180 and AFGL 5142. 



-29- 



(a)MonR2 a 



Si *■ 



m^'. :■! 



B -34' 



Al>Ofl'"ID' 



fttfNW* ai™4Q» or^jC 

Riglit AatieiiKion (J£ODD> 



*v*"'ao* 




(b) Serpens 



^■;.^^ 









AnjEHlnr Of^TvL 




-6 D' 



^;t-^ 



aij«":c aiif"QO" a'^hu' 2fl"«« 



Fig. 10. — Same as Figure |8] for Mon R2 (a) and Serpens SVS2 (b). The contours in the 
P-V diagrams with the intervals of the 3 a levels start from the 3 a levels, where the 1 a 
noise levels are 0.07 K, and 0.08 K in T^ for Mon R2 and Serpens SVS2. 
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Fig. 11. — The maps of blue-shifted and red-shifted velocity range of the H^^C0^(J=1- 
0) emission (contours) superposed on the JHKs composite color images in log scale from 
2MASS for S87E (a), S88B (b), AFGL 5142 (c), and AFGL 5180 (d). The contours with 
the intervals of the 3 a levels start from the 3 a levels. The 1 a noise levels for the blue- 
shifted range are 0.09Kkms~^, 0.06Kkms~^, 0.05Kkms~^, and 0.04Kkms~^, and those 
for red-shifted range are 0.07Kkms~^, 0.06Kkms~^, 0.08Kkms~^, and 0.07Kkms~^in T^ 
for S87E, S88B, AFGL 5142, and AFGL 5180, respectively. The white arrows show 1 pc 
scale of individual regions. The velocity ranges of blue-shifted components are 20kms^^ 
to 22kms^"'^, 20kms^"'^ to 22kms^"'^, — 5kms^"'^ to — 3kms~^, lkms~^ to 4.2 km s""*^ for 
S87E, S88B, AFGL 5142, and AFGL 5180. The velocity ranges of red-shifted components 
are 22kms^^ to 27kms~^, 22kms~^ to 23kms~^, — lkms~^ to — 3kms~^, 4.2kms~^ to 
5.4kms-i for S87E, S88B, AFGL 5142, and AFGL 5180. 



Table 1. Physical parameters of the H CO clumps 



Source Name 


RA (J2000) 


Dec (J2000) 


D[pc] 


Mcluster [Mq\ 


Kclump[pc] 


AVciumpIkm s-i] 


Af(H2)[x 1022 cm-2] 


Mclump[MQ] 


Tex [K] 


S87E 


19:46:19.9 


24:35:24 


2100 


180 


0.75±0.4 


4.0±0.4 


4.0±0.2 


1400±400 


22 


S88B 


19:46:47.0 


25:12:43 


2000 


120 


0.37±0.2 


2.7±0.4 


4.8±0.3 


420±120 


32 


Gem 4 


06:08:41.0 


21:30:49 


1500 


190 


0.58±0.3 


2.0±0.3 


4.5±0.2 


940±200 


28 


AFGL 5180 


06:08:54.1 


21:38:24 


1500 


60 


0.43±0.2 


2.4±0.3 


4.3±0.2 


480±100 


24 


GGD 12-15 


06:10:50.9 


-06:11:54 


830 


73 


0.44±0.1 


2.0±0.3 


2.9±0.2 


340±50 


19 


S235AB 


05:40:52.5 


35:41:25 


1800 


220 


0.61±0.3 


1.7±0.3 


2.3±0.1 


520±130 


24 


AFGL 5142 


05:30:45.6 


33:47:51 


1800 


50 


0.54±0.3 


3.6±0.3 


3.1±0.2 


550±140 


17 


Mon R2 


06:07:46.6 


-06:22:59 


830 


340 


0.75±0.2 


1.9±0.3 


2.8±0.2 


930±140 


19 


S140 


22:19:18 


63:18:48 


900 


12 


0.42±0.1 


2.0±0.3 


6.3±0.4 


690±100 


24 


AFGL 490 


03:27:38.7 


58:46:58 


910 


25 


0.45±0.1 


1.8±0.3 


2.1±0.1 


250±40 


16 


BD-|-40°4124 


18:29:56.8 


01:14:46 


900 


12 


0.30±0.1 


1.5±0.1 


1.9±0.1 


100±20 


15 '^ 


NGC 7129 


21:43:02 


66:06:29 


1000 


76 


0.48±0.2 


1.5±0.4 


1.4±0.1 


140±30 


12 


Serpens SVS2 


18:29:56.8 


01:14:46 


260 


27 


0.24±0.1 


1.7±0.3 


3.0±0.2 


llOitlO 


14 





Note. — This table shows the name, center coordinates, distance, mass of the clusters (from Lada & Lada (2003), Porras et aL (2003)), Rdump- radius, AVdump- 
velocity width, N{R2)'- mean H2 column density, Mdump: clump mass, Tox: excitation temperature derived by kinetic temperature of NH3 of the clumps. 

^Because wc could not estimate the excitation temperature for low signal-to-noise ratio of the NH3 data, we adopted 15 K which is the typical temperature in 
intermediate regions. 
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Table 2. Score sheet of the morphological classification of H CO clumps 



Source Name 


Type 


SFE 


The Number of 


Brightness 


Offset Between 


Aj 


^Kclump 


(H13CO+) 


Md,n,p(Hl3CO+) 






[%] 


Sub-clumps 


Contrast 


Clumps and Clusters [pc] 


[pc] 


/Ay,i,^p{ciso) 


/M,i,^p(Ci80) 


S87E 


B 


12±3 


2 


High 


Off-peak : 0.3-0.6 


0.17 


1.2±0.1 




0.85±0.05 


S88B 


C 


22±5 


1 


Low 


Off-peak ; 0.5 & Cavity : ~ 0.5 


0.13 


0.8±0.1 




0.09±0.01 


Gem 4 


B 


17±3 


2 


High 


Peak & Off-peak : 0.9 


0.16 


l.liO.l 




0.75±0.04 


AFGL 5180 


B 


11±2 


1 


High 


Off-peak : 0.2 


0.13 


l.OiO.l 




1.0±0.26 


GGD 12-15 


B 


18±3 


> 2 


Low 


Off-peak ; 0.3-0.5 


0.14 


0.9±0.1 




0.25±0.01 


S235AB 


C 


30±6 


> 2 


Intermediate 


Off-peak : 0.3-0.5 


0.21 


0.8±0.1 




0.40±0.07 


AFGL 5142 


A 


8±2 


1 


High 


Peak 


0.14 


1.9±0.1 




0.58±0.03 


Mori R2 


C 


27±3 


> 5 


Low 


Off-peak : 0.2-0.5 & Cavity : ~ 0.3 


0.19 


0.9±0.1 




0.77±0.07 


8140 


A 


3±1 


1 


High 


Peak 


0.11 


0.8±0.1 




0.37±0.01 


AFGL 490 


B 


9±1 


3 


Intermediate 


Peak & Off-peak : 1-1.2 


0.16 


0.8±0.1 




0.13±0.01 


BD+40°4124 


B 


10±2 


2 


Low 


Peak & Off-peak : 0.1-1 


0.13 


1.1±0.2 




0.49±0.03 


NGC 7129 


C 


35±5 


> 4 


Low 


Off-peak ; 0.2-1.5 & Cavity : ~ 0.5 


0.20 


1.0±0.2 




0.19±0.01 


Serpens SVS2 


C 


20±1 


> 4 


High 


Off-peak : 0.2-0.5 


0.09 


0.9±0.2 




0.71±0.04 



00 
to 



Note. — This table shows the source name, the classified type of the clump, the SFE derived by using the H^'^CO^ mass, the number of sub-clumps, the 
brightness contrast of the clump, the location of the cluster, the Jeans length, the velocity width ratios of H^'^CO"'" to C^*0, and the mass ratios of H^^CO"'" to 
C^*0. The location of the cluster is written as follows. Peak: The emission peak of the clump which is associated with the cluster; Off-Peak: The emission peak of 
the clump which is not associated with the cluster. The offset scale between the peak and the cluster is subsequently-shown [pc]; Cavity; In case that the cluster 
is located in a cavity. The size of the cavity is subsequently-shown [pc] . 



Table 3. The virial parameters of the H CO clumps 



Source Name 


/3rot 


Aurb 


Ahcrm 

[xlO-*^] 


/3vir 


O-grad 

[km s-1 pc-1] 


-^•^star — max 

[Mq] 


References 


S87E 


0.79±0.25 


0.53±0.16 


2.4 


1.6±0.51 


4.3±0.9 


20±4 


1 


S88B 


0.34±0.11 


0.56±0.17 


5.8 


0.81±0.25 


2.2±0.4 


20±4 


2 


Gem 4 


0.05±0.01 


0.23±0.06 


3.5 


_0.44±0.10 


1.0±0.2 


19±4 


3 


AFGL 5180 


0.33±0.08 


0.47±0.11 


4.4 


0.6±0.14 


2.2±0.4 


18±4 


4 


GOD 12-15 


0.05±0.01 


0.48±0.08 


4.9 


0.07±0.01 


0.7±0.1 


15±3 


5 


S235AB 


0.09±0.02 


0.29±0.08 


5.7 


-0.24±0.07 


1.0±0.2 


15±3 


6 


AFGL 5142 


0.55±0.15 


0.91±0.25 


3.2 


1.9±0.53 


2.7±0.5 


10±2 


7 


Mon R2 


0.18±0.03 


0.13±0.02 


3.1 


-0.38±0.07 


1.7±0.3 


10±2 


8 


S140 


0.05±0.01 


0.22±0.04 


3.0 


-0.46±0.08 


1.0±0.2 


10±2 


9 


AFGL 490 


0.20±0.03 


0.28±0.05 


5.9 


-0.04±0.01 


1.2±0.2 


8±2 


10 


BD+40°4124 


0.05±0.01 


0.55±0.11 


8.7 


0.22±0.04 


0.5±0.1 


8±2 


11 


NGC 7129 


0.23±0.05 


0.55±0.12 


8.5 


0.56±0.13 


0.9±0.2 


8±2 


12 


Serpens SVS2 


0.33±0.03 


0.60±0.05 


6.5 


0.85±0.08 


1.4±0.3 


3±1 


13 





00 
00 



References. — (1) Chen ct al. (2004) ; (2) Garay et al. (1993) ; (3) Ghosh et ah (2000) ; (4) Minier et ah (2005) ; (5) 
Fang & Yao (2004) ; (6) Felli et ah (1997) ; (7) Zhang et ah (2002) ; (8) Carpenter (2000) ; (9) Preibisch & Smith(2002) 
; (10) Schrcycr ct ah (2002) ; (11) Looncy ct ah(2006) ; (12) Fucntc ct ah (2002) ; (13) Kaas ct ah (2004) 



Note. — This table shows that the the virial ratios of rotation, turbulent, thermal, and the combination of the above 
energies, which are derived adopting p = 0, respectively, the velocity gradients within the clumps, the highest stellar 
mass of the cluster members, and the references of the highest stellar mass. 

^The highest stellar mass of the cluster members. 



